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Abstract 

We analyze the production of </> mesons in vrA^ and A^A^ reactions in the near- 
threshold region, using throughout the conventional "non-strange" dynamics 
based on such processes which are allowed by the non-ideal to — (/) mixing. 
We show that the occurrence of the direct (pNN interaction may show up 
in different unpolarized and polarization observables in vrA^ N(p reactions. 
We find a strong non-trivial difference between observables in the reactions 
pp — > pp4> and pn pn<j) caused by the different role of the spin singlet 
and triplet states in the entrance channel. A series of predictions for the 
experimental study of this effect is presented. 
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I. INTRODUCTION 



The present interest in the meson production in different elementary reactions is related 
to the strangeness degrees of freedom in the nucleon. Since the meson is thought to consist 
mainly of strange quarks, i.e. ss, with a rather small contribution of the light u and d quarks, 
its production should be suppressed if the entrance channel does not possess a considerable 
admixture of strangeness. Indeed, the recent experiments on the proton annihilation at rest 
(cf. for references and a compilation of data) point to a large apparent violation of the 
OZI rule [0, which is interpreted as a hint to an intrinsic ss component in the proton. 
However, the data can be explained as well by modified meson exchange models |^ without 
introducing any strangeness component in the nucleon or OZI violation mechanisms. On the 
other hand, the analysis of the vrA^ sigma term suggests that the proton might contain a 
strange quark admixture as large as 20%. Thus this issue remains controversial. Therefore it 
is tempting to look for other observables that are sensitive to the strangeness content 

of the nucleon. Most of them are related to a possible strong interference of delicate ss 
knock-out and shake-off amplitudes and the "non-strange" amplitude which is caused by 
OZI rule allowed processes, or by processes wherein the standard OZI rule violation comes 
from the cj) — uj mixing. 

As shown in Ref. 0, through this interference the polarization observables of the 
photoproduction process are sensitive even to a rather small strangeness admixture in the 
proton. However, the only ^'^S ss configurations may be seen in a such process. The 
other configurations, such as ^'^P, are suppressed by the selection rules and/or form factors. 
Contrary to this, Ellis et al. argue that the possibly dominant ^Pi configuration might 
be seen in hadronic reactions. 

Obviously, reliable information about the hidden strangeness manifestation in the 
production in irp and A^A^ reactions can be obtained only when the conventional, i.e. non- 
exotic, amplitude has been understood quantitatively. This is the objective of this work. 
The dominant conventional processes in ttN and A^A^ reactions are depicted in Figs. 1 and 2, 
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where (a) is the mesonic exchange process being allowed by the finite (ppir coupling strength 
and (b) is the direct (f) radiation from the nucleon legs, which is proportional to the finite 
(f)NN interaction strength. It should be emphasized that the process la is a subprocess in 
the diagram 2a, while the process lb is a subprocess in the diagram 2b, when the exchanged 
boson is a pion. 

While the diagrams in Figs. 1 and 2 look like usual Feynman diagrams it should be 
stressed they give a guidance of how to obtain from an interaction Lagrangian of hadronic 
fields a covariant parameterization of observable in strict tree level approximation. Ad- 
ditional ingredients are needed to achieve an accurate description of data within such a 
framework. In particular, the vertices needs to be dressed by form factors. 

The early theoretical studies show, indeed, that predictions for hadronic observables 
are very sensitive to the parameters of the monopole form factors which can not be fixed 
unambiguously without adjustments relying on the corresponding experimental data. In 
our case one can rely on the recent measurement of the ratio of the total cross sections of 
and uj production in pp reactions studied by the DISTO Collaboration at SATURNE at 
Tiab = 2.85 GeV as well as on the angular distribution in the pp ppcf) reaction and 
on the total cross section |T] 



In the (sub)threshold region also in heavy-ion reactions the production data is accessible 
via the K^K~ decay channel studied with the 47r detector FOPI at SIS in GSI/Darmstadt 



12 1 . However, here an understanding of elementary hadronic reactions serves as a prereq- 



uisite for interpreting the data. Rather the upcoming proton and pion beam experiments 



with the HADES detector system at SIS in GSI |jT3|] offer a chance to enlarge the data base. 
In particular, HADES can identify the also via the e^e~ channel. 

Finally we mention that for pion-induced reactions at the proton also near-threshold 



data for the production are available |]T^| 



An important step towards an understanding of the structure of the "non-strange" 
pp —* pp(j) reaction mechanisms was made recently . The focus of Ref . |]r3| is the determi- 



nation of the parameters of the direct (pNN interaction, thus reducing the above mentioned 



uncertainties, by analyzing the combined the pp ppu and pp ppcj) reactions and the 
corresponding DISTO data ||I0| at a given beam energy; for the reaction pp —>■ ppuj just the 
same mechanism is assumed as those for pp —>■ pp(j), as shown in Fig. 2. Clearly, at one given 
beam energies the excess energy for both reactions are quite different. 

In this paper we therefore attempt a different approach with a similar goal. We reduce 
the uncertainties of the reaction mechanism by a combined study of the to each other 
related reactions pp ppip and Ti^p nip using the known data within the same interval 
of excess energies of 20 — 100 MeV [0,0 • For the absolute normalization of the angular 
distribution in pp ppcj) ||10| we use the recently published total cross sections of the reaction 



pp — > ppcj) |Tl|]. It turns out that this value is compatible, within given error bars, with an 
extrapolation of the previously measured ratio of the total cross sections of pp — > ppcj) to 
pp ppcj) reactions, a^/a^^ = (3.7±0.5) x 10^^ at Tiab = 2.85 GeV ||10|, by normalizing it to 



the old bubble chamber data for at various excess energies [16|. The extrapolated value 



of the cross section differs from the new value |Tl[] thus influencing to some extent adjusted 
parameters. 

In comparison with previous works we are doing the next step towards an understanding 
of the dynamics of cj> production in hadronic reactions. We present a systematical analysis 
of TT^p — * ncj), pp ppcj) and pn —* pncj> reactions in the near-threshold region where 
the destructive interferences between the two mechanisms (a) and (b) in Figs. 1 and 2 are 
essential. We are going to show that basically there are two different sets of the model 
parameters. One of them corresponds to the case when the mesonic exchange channel (a) 
is dominant for the n^p — > ncj> reaction (Fig. 1), and in the second case the direct emission 
mechanism (b) is dominant. For both sets of solutions we calculate the total and differential 
cross sections and spin density matrix, responsible for the cj> e+e^ and cj> K^K^ decay 
angular distributions and show in which observables the direct cpNN interaction might be 
clearly manifest. We present also a combined analysis of pp ppcj> and pn pncj> reactions 
at a finite energy excess with taking into account the final state interaction and analyze the 
deviation of predicted observables from the pure threshold values which is important for the 



future understanding of the role of an internal strangeness degrees of freedom in the nucleon. 
For this aim we study the beam-target spin asymmetry and the relative role of the singlet 
and triplet states in the entrance channel. 

Our paper is organized as follows. In Section II, we define the kinematical variables 
and formulae for calculating the cross sections and polarization observables. The basic 
amplitudes for the mechanisms illustrated in Figs. 1 and 2 are given explicitly in Section III. 
In Section IV we discuss results of numerical calculations and predictions. The summary 
is given in Section V. In the Appendix we describe the formalism for the enhancement 
factors of final state interaction within the framework of the Jost function and the effective 
phase-equivalent potentials. 



II. OBSERVABLES 

The differential cross section of the reaction 7r~p — > 4>n (cf. Fig. 1) has the obvious form 
where p^^ = {E^^, p^) and q = {E^, q) are the four-momenta of the pion and the 4> meson in 



the center of mass system (c.m.s.); means average and sum over the initial and final 

spin states, respectively. 

The differential cross section of 4> production in the reaction a + b ^ c + d + 4>, where 
a, b and c, d label the incoming and outgoing nucleons (cf. Fig. 2), is related to the invariant 
amplitude T(2) as 

da = / % % ^ Si'^ (P. -Pf). (2) 

2(27r)V«(«-4M2) ' ^ 2E,2Ed2E^ \ ^ i) \ J 

where p„ = {E^, Pn) with n — a,b,c,d are the four-momenta of the nucleons in the c.m.s., 
^/s = Ea + Ef) is the total c.m.s. energy, Pjj are the total four-momenta of the initial or 
final states. Hereafter 9 denotes the polar meson angle and Q is its solid angle. We 
use a coordinate system with z || p^, y || PaXq. Among the five independent variables 
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for describing the final state we clioose E^, Q and flc- Then the energy Ec of particle c 
is defined by E, = ^^-gV^^-g(^! z^^ ^-^^ ^ ^ 2(v/i - E^), B = s- 2E^^s + M|, 
C = 2|q| cos6'qp^, Finally, the fivefold differential cross section reads 
c/V 1 





|q| 


IpcI 


2 


\A\ 


Pel 


+ CE,\ 



^(^) ^.iLi'"^^.i - (3) 



dE^dVtdVt^ %{2Txf ^ s{s - AMI) 
The total and/or partially differential cross sections are found by integration over the avail- 
able phase space. 

In this paper we consider two polarization observables. One of them is the spin density 
matrix which describes the spin structure of the outgoing meson, 

Erji rji^ 

Prr' sr^ rji rji^ 1 \^) 
2-^r,l3 ^r,l3 r,l3 

where r = = ±1,0 are the spin projections of the meson, and (3 denotes a set of 
unobserved quantum numbers. The spin density defines the angular distribution in — »• 
e^e~ and (f) —>■ K^K~ decays, which has a simple form in a system where the (f) meson is at 
rest (for details see [^). The decay angles O, $ are defined as polar and azimuthal angles of 
the direction of flight of one of the decay particles in the meson's rest frame. The decay 
distributions integrated over the azimuthal angle $, VV(cos6), depend only on the diagonal 
matrix elements poo; Pii = P-i-i? normalized as poo + 2pii = 1, according to 

>V(cos 6) = ^(^^ {^ + B cos' e) , (5) 

where the decay anisotropics B read 

j^K+K- ^ 1 ~ 3P00 j^e+e- ^ ^ ~ ^Pqq 

1 - Poo ' 1 + Poo ' 

To exclude the kinematical dependence of poo or B on the meson production angle, we 
choose the quantization axis along the z direction (in the rest system), and using the 
corresponding Wigner rotation functions d'^^.ix) o^is gets the amplitudes T^^p in Eq. @) by 



where only = —9 is important, while the other x^'s disappear in Eq. (^. 

Another polarization observable is the beam-target asymmetry in the NN NNcf) 
reactions which is related to the nucleon spin states via 

^ dajSj = 1) - da{Si = 0) 

~ da{Si = 1) + daiSi = 0) ' ^ 

where Si is the total spin in the entrance channel. It is important to note that spin and 
parity conservation arguments result in a precise model independent prediction fl^ for the 



beam - target asymmetry: Cbt = ^ for the pp ppcf) reaction at the threshold. In the 
pn — >■ pncf) reaction the asymmetry depends on the relative weights of the triplet and singlet 
states in the entrance channel. 

III. BASIC AMPLITUDES 

Basically, our consideration in this section is similar to the previous study p| (for the 



pure mesonic exchange contributions depicted in Figs, la and 2a) and to the models PJT5 
for both channels shown in Figs. 1 and 2. The difference between this work and previous 
ones is in the different form of cut-off form factors for the off-shell nucleons in direct 
emission (cf. Figs, lb and 2b) and a different choice of the cut-off parameters in vriV and 
A^A^ interactions which we will discuss below in detail. In spite of the mentioned similarity, 
for completeness in discussing our predictions for the set of observables which have not 
been considered before, in this section we display the main formulae which define the basic 
amplitudes. The meson - nucleon and the 0p7r interaction Lagrangians read in standard 
notation 

CmNN = -ignNNN-f^TTvN 

C^p. = g^pn e^""'' 9^0. Tr(9„p^7r), (10) 
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where Tr(p7r) = p^vr" + p+vr + p vr"*", and bold face letters denote isovectors. All cou- 



pling constants with off-shell meson are dressed by monopole form factors |[T^ Fi = 
(Af — 'mf)/{Aj — kf), where fcj is the four-momentum of the exchanged meson. Follow- 
ing the scheme of the meson photoproduction [jl9| we assume that (pNN vertices must be 
dressed by form factors for off-shell virtual nucleons. But this might result in a violation of 
the transversality of the amplitude with respect to the meson field. To avoid this problem 
we use the prescription of Ref. |jl9| and parameterize the product of the two form factors 



appearing in the left and the right diagrams in Figs, lb and 2b in a symmetrical form 

1 f A% A% \ 

- 2 [a% + ip^ - M,Y + a% + [vI-m,y) ' 

here pi (pr) is the four-momentum of the virtual nucleon in the left (right) diagrams in 
Figs, lb and 2b. 



A. 7vN — > Ncj) reaction 

The invariant amplitude for the meson exchange channel (a) in Fig. 1 reads 

T(ia)A = K-^'e'^'^ [u{p,)rpi{k,)u{pa)] q^kkefI^, (12) 

where 

r;(A:) = y + ^^a'^A:,,, (13) 

A 2 ^2 AP2 _ 2 

T^-^N(j^ \ _ 9pNN g4>fm ^^pNN ^^<t>pn "^p 

S ^p ^^pNN % A^^^ - 

with kp = pc — Pa as the virtual p meson's four-momentum; is the meson polarization 
vector, denotes the isospin factor being equal to (1) for a vr^ (n^) meson in the 
entrance channel, and the nucleon spin indices are not displayed; i,j--- are Lorentz indices, 
and 7j and u denote Dirac matrices and bispinors. 

The invariant amplitude for the direct radiation channel (b) in Fig. 1 has the following 
form 

8 



T{lb)\ = WcjiNN 9ttNNu{Pc 



u{pc)e*^ I^FNipi,Pr), (15) 



where r^(g) and are defined by Eqs. (|13]) and (|TT]), respectively, and pi = Pa — Q and 
Pr=Pc- q- 

B. AT AT NNcf) reaction 

The total invariant amplitude of meson exchange diagrams (a) in Fig. 2 with internal 
meson conversion is the sum of 4 amplitudes 

(TA/)a = CTm [ab; cd] + ^Tm [ab; dc] + ^Tm [ba; dc] + ^^Tm [ba- cd] (16) 

with = = -ilp = -^^p = 1, = = -1, = = -2. The last two terms 
stem from the antisymmetrization or from charged meson exchange in pp or pn reactions, 
respectively^ The first term in Eq. (|16[) for the pp reaction reads 

^2a)x[ab- cd\ = K^^ [u{p,) 75 u{pk) ] [u{p,) Ti{k)uipa) u,Mk; qiel'] , (17) 

with 

T^NNf,2 , 2n ^ 9nNN QpNN g<t>piT - Ap - A^^^ - m?p A^^^ - ml 

^ ^' (kl - mlWp - ml) Al -kl k^- Af^^ - k^ Kf^^ -k^' ^ ' 

The amplitude of direct meson emission from the nucleon legs according to Fig. 2b is 

calculated similarly to the real or virtual photon bremsstrahlung |2^j2l| in the few GeV 

region. The internal zig-zag lines in Fig. 2b correspond to a suitably parameterized A^A^ 

interaction in terms of an effective two-body T-matrix which is written in the form of the 

one-boson exchange model (OBE) with effective coupling constants and cut-off parameters 

and may be interpreted as effective vr, u, p, a meson exchanges. We would like to stress 



^In p we used a convention with = = 2, which, however, does not change our threshold 
prediction for the ratio of the total cross sections in pn and pp interaction made there without the 
final state interactions. 



that this is an effective dynamical model which is appropriate in the few GeV region and 
which is different from the OBE model in the conventional sense. This model has been 
applied successfully to different reactions pO|-p^ and this encourages us to employ it for the 



(j) production too. 

The total amplitude for the process (b) in Fig. 2 consists of 2-8 (2-6) contributions for pp 
(pn) interactions and has a similar structure as Eq. (|16D (with = = 1, = = 0, 
for (T,u! exchanges), where T[ab; cd] now reads 



T(2b) \[ab; cd] = -g^NN^*^ [u{pd)V'^u{pb)] x 



J2 Hd^]Hpc 



^^F^ r;(-g) + r;(-g) ^ v 



u{Pa). (19) 



m=7r,cr,p,a; 

Here, and D*" are effective coupling vertices and propagators of the two-body T matrix, 
respectively, 

where k is the four momentum of the virtual meson m and G^nn is the vertex function 
which includes the corresponding cut-off form factor. The numerical values of the GmNN 
are taken from Refs. |20|j2^ . 

In the near-threshold region the relative velocity of the outgoing nucleons is small which 
might result in a strong final state interaction (FSI) between them. If the energy excess is 
a few MeV up to a few tens MeV then one can consider only the s-wave interaction and 
account for the final state interaction in terms of the enhancement factors by renormalizing 
the basic amplitude correspondingly. For instance, for the pn reaction we get 

) , (22) 

where mcmd are the spin projections of the nucleons in the final state, and Xq, Xi are the 
singlet and triplet enhancement FSI factors, which are calculated within the Jost function 
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and the phase-equivalent potentials formalism, which we describe in detail in Appendix A. 
The calculation shows that the singlet enhancement factor is much greater than the triplet 
one, i.e. |Xop — 1 ^ |Xip — 1 ^ and greater than the corresponding factors of higher 
partial waves. Thus, for the pp interaction we use 

Tpp[ab] cd] Tpp[ab; cd] (Jopp 6 

) , (23) 

reminding that at the threshold the pp triplet final state is exactly zero, lopn and lopp '^^^ 
different which reflects the difference in the corresponding effective radii and the scattering 
lengths. Note that the mutual FSI of the (p and the outgoing nucleons is assumed to be 
negligible. 

In calculating the cross sections and the spin density matrix, squares and bilinear forms 
of the FSI-corrected amplitudes need to be evaluated numerically. 



IV. RESULTS 



A. Fixing parameters 



The parameters of the two-body T matrix for the direct emission depicted in Fig. 2b 
are taken from Refs. P0J22], where a quite reasonable agreement with data of different elastic 



and inelastic A^A^ reactions is found. 

The coupling constant g^pp^, is determined by the (j) ^ pn decay. The recent value 
r(0 pn) = 0.69 MeV results in \g(f,pn\ = 1-10 GeV^^. The SU(3) symmetry consideration 
T5| , |23[| predicts a negative value for it. Thus, g^ppn = —1.10 GeV^^. 

The remaining parameters of the meson exchange amplitudes for the processes in Figs, la 



and lb are taken from the Bonn model as listed in Table B.l (Model II) of Ref. |T8 

The yet undetermined parameters are: the cut-off parameters for the virtual mesons in 
the (ppn vertex, AJ^ and A^^, the cut-off A^r in Eq. (0), and the parameters of the (pNN 
interaction, g^/^NN and k^. We can reduce the number of parameters by making the natural 
assumption AJ^ = A^^ based on the symmetry of the virtual mesons in the (ppir vertex []T3 
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The next consideration is related to the tensor couphng k,^. Based on the cf) — uo similarity 
we do not expect a large value for it and in all our subsequent calculations we employ the 



theoretical estimate = 0.2 as an upper limit. 

Even after that we have three free parameters being gt/^NN, and A^r. For g^NN the 
SU(3) symmetry predicts p5[ ] 

g^NN = -tg'^Ov gujNN, (24) 

where AOy is the deviation from the ideal u — (p mixing angle. It is responsible for the 
"standard" OZI rule violation, and in general, it depends on the method of its determination 
(Gell-Mann-Okubo linear or quadratic mass formulae, radiative decays, say (f){uj) 771, 
etc.). Using the quadratic Gell-Mann-Okubo mass formula one gets AOy = 3.7°. Sometimes 
another relation is used, e.g. 

g<t>NN = -3 sin A9v gpNN, (25) 
which is obtained from the SU(3) relation 

guiNN = jp^jj COS A9v gpNN, (26) 

and Eq. ([2^) with the assumption D/F = in the SU(3) vector meson octet. Using the 
known values for glj^^/'^'^ — ^-"^ ~ ^-^ HH and gl^N^/'^T^ = 22 — 24 |T^, one may obtain 



—g(j>NN = 0.57 — 0.65 and —g^NN = 1.07 — 1.09 for the expressions ( p5D and p^), respectively. 
On the other hand, the theoretical estimates of Ref. give g(j>NN = —0.24. Thus, we can 



conclude that even using the standard OZI rule violation (thought non-ideal u — (p mixing) 
one is left with estimated values of g^^NN within a quite large interval. The possible hidden 
strangeness in a nucleon may even increase this interval. In this paper we restrict ourselves 
to the standard OZI rule violation mechanisms and analyze consequences of varying —g^NN 
in the region 0.0 - 1.0. 

The negative coupling constant grppj^ results in a destructive interference between meson 
exchange amplitudes (a) and direct emission (b) in Figs. 1 and 2. Analyzing the unpolarized 
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7T p ^ ncj) reaction, based on the data of Ref. |]14|, we find that the yet unconstrained three 



parameters q^nni A^^, A^v become related to each other as ^1p=^1p{,g(j,NN^ A^^) by the 
constrains given by the data, and two solutions emerge for this dependence: (i) cr(a) > cr(b) 
and (ii) o"(a) < cr(fe), where o't^a,h) are the total cross section for the meson exchange process 
(a) and the direct emission (b) calculated separately. These solutions are displayed in 
Fig. H for several values of q^nn as discussed above. 

In order to constrain one more free parameter we analyze also the cross section da/dfl 



for the pp ppcf) reaction, using simultaneously the DISTO data [I^JTll, i.e. the angular 
distribution [|l^ (we remind that in our notation VL is the (f) meson solid angle) and the total 
cross section [|ll| , respectively. For this aim we fix the absolute normalization of the angular 



distribution da/dQ given in Ref. ||T0[ by making use of the recently published the total cross 



section [TT[]. As a result we get the fat dots in Fig. ^. The minimum values for —g^f^NN are 



0.07 (for > <7(^b)) and 0.60 (for (T(a) < <J(b)), respectively at — > oo. For both solutions 
we find that increasing values of \g(f,NN\ results in decreasing values of A^r leaving the total 
cross section or "integrated" strength of the (pNN interaction on the same level. 

To explore in more detail the relative importance of the direct (pNN interaction we now 
employ three parameter sets; two of them correspond to the (T(a) > (7(6) solution and the 
third one to a^a) < <7{b)- 

set A: g^N^ = -0.24 [H, A^^ = 1.34 GeV, A^ = 1.065 GeV (a(„) > a^i,)), 



set B: g^NN = -0.8, A^^ = 1.34 GeV, A^ = 0.715 GeV (a^^) > (^(1,)), 
set C: g^NN = -0.8, A^^ = 1.60 GeV, A^ = 1.99 GeV (a(,) < a^b))- 

In the sets B and C we choose —g^jjNN close to its upper limit as predicted by SU(3) symmetry. 
Results of our calculation of the total cross section for the 7r~p — ^ n0 reaction for these 

1 1 i 

parameter sets are shown in Fig. ^ as a function of the energy excess As2 = — Sq^j^ 
with sottat = {M]y + m^Y. The contribution of the meson exchange channel (a) is displayed 
separately by the dot-dashed line, and the direct emission by the nucleon line (b) is depicted 
by the dashed line. Clearly seen is the strong destructive interference of the channels (a) 
and (b), in particular for the set C. 
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In Figs. |] and ^ we display results of our calculations for da/dfl for the pp ppip reaction 



together with the available data |jTO|,|TT| for the parameter sets B and C, respectively, at 

1 1 i 

As2 = S2 — sqIj^ = 82 MeV with Sqatat = (2M7V + m<^)^. The results for the set A are very 
similar to that of the set B and we do not separately display them here. Interestingly, in all 
parameter sets considered the channel (a) dominates in the NN —>■ NNcp reaction, but the 
interferences are different for different parameter sets. 

One can see a qualitative difference in irN and A^A^ reactions for the set C. In the vrA^ 
reaction the direct radiation channel (b) is dominant, i.e. cr(^h) > <^{a)- In contrast, in the A^A^ 
reaction the relative contribution of the direct emission channel (b) increases considerably as 
compared with sets A and B, but it is still smaller that the meson exchange channel (a). The 
reason for this difference is the following one. The direct emission (b) in the reaction vrA^ 
consists of the two competing u-channel and s-channel diagrams shown in Fig. lb. which add 
destructively, while the contribution of the of the u-channel amplitude is greater. However, 
the corresponding contributions of the two competing diagrams in Fig. 2b are numerically 
nearly the same resulting in a stronger suppression of the direct channel (b). 

As we have adjusted our parameters by the data, it is clear that they describe the data 
with approximately equal quality, and at the present level of the data accuracy it is difficult 
to give a preference to one of them. Therefore, we now investigate whether other observables 
can be used to constrain the parameters further and whether the difference between pp and 
pn reactions is a sensible measure. 



B. TT p — > n<j) reaction 

The calculated angular differential cross sections of the 7r~p ncj) reaction at As^ = 50 
MeV and for the parameter sets A, B, C are shown in Fig. ^ One can see that the shapes of 
the distributions for the sets A and B are very similar to each other. They are quite smooth 
and close to the distribution of the meson exchange channel (a). Only in the backward 
direction the total cross section slightly decreases due to the destructive interference with 
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the direct channel (b), leading to some enhancement of the cross section in forward direction. 
(We would like to mention here, that the extrapolated value of the (f) production cross section 



in A^A^ reactions |jT6[| from the data [|T0| would require somewhat different parameter sets, 
which in turn cause also more pronounced differences between the sets A and B.) Contrary 
to that, for the model C the largest destructive interference appears at forward direction, 
where the contributions of the two competing channels (a) and (b) are close to each other. 
As a result, the cross section is enhanced in the backward direction. So we can conclude 
that the differential cross section is sensitive to the dynamics of the production and the 
direct (pNN coupling (a similar conclusion for uj production has arrived at in Ref. [p^] ). 



The prediction for spin density matrix elements poo Pii for the different models is 
shown in Fig. ^ as a function of the production angle in c.m.s. at As^ = 50 MeV. The sets 
A and B deliver standard values, typical for the spin-flip processes, i.e. poo — 0, pn — 0.5. 
But the parameter set C predicts a strong deviation from these values, especially at forward 
direction. The reason of this effect is the following. In the meson exchange channel (a) 
the nucleon spin-flip amplitudes result in transitions — m^,mj, where ruf = m.i — 
with = ±1. For instance, the transitions like — | —1, | are dominant. In the direct 
radiation channel (b) together with this strong amplitudes we have finite amplitudes for the 
transition rrii — 0, mj, where fnf = rrii. In the set C the strongly competing nucleon spin flip 
amplitudes cancel each other and only the nucleon spin conserving direct emission amplitude 
(b) survives. This is illustrated in Fig. ^, where we show the nucleon spin flip (left panel) 
and the nucleon spin conserving (right panel) amplitudes for the set C. Here = ImT'^P^'"'^ 
with the quantization axis along the z direction (T'^p^^'I' is purely imaginary). 

The anisotropics of the decay channels e~^e~ and K^K~ (cf. Eqs. (|, P) for 



the different parameter sets are shown in Fig. 10. Again, one can see a strong deviation 



of our prediction for the set C from the naive expectation 5'^^'^ ~ 1, B^^^ ~ —1 based 



on a purely mesonic exchange channel or on the sets A and B. Fig. ^ illustrates the the 
manifestation of this deviation in the real e~^e~ and K^K~ angular distributions. The 
distributions Wl and Wt are the longitudinal (along the quantization axis) and transversal 
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fluxes for tlie outgoing electrons or kaons. The functions Wl and Wt are normalized as 
/ a/M/|(0) + W^{Q)dcosQ = 1, where 6 is defined by Eq. (|]). Thus, one can see that the 
sets A and B predict a practically vanishing kaon flux in the longitudinal direction for all (j) 
production angles. The set C predicts a flnite amount of the longitudinal flux which increases 
with decreasing (p production angle in c.m.s. A corresponding modiflcation is predicted for 
the electron flux too. 



C. pp — > Pp4> ^nd pn — > pncj) reactions 

As shown in Figs. ^ and ^ the meson exchange contribution (a) is the dominating contri- 
bution to the NN — s> NN(f) reaction, therefore, it is useful to recall the threshold prediction 
for this channel in the absence of the flnal state interaction, which serves as a starting point 



for further calculations at flnite energy. Adopting the notation of Ref. [0 we can express 



the invariant amplitudes of the reaction ab ab(f) with a = p, b = p or n as following 

Tpp = Fi, Tpn = - (Fq + Fi) , (27) 

where Fq (Fi) is the initial singlet (triplet) amplitude with 

i^O = /o(-l)^+™" S-mamA^lm. ^Im, ' <^-imJ, 

Fi = /i(-1)H™^ 5™.„.,(5|^^. - 5i^J, (28) 

where rria^b and m^^ are again the nucleon spin projections in the initial and flnal states, 
respectively, /o = Gy^To, /i = 2\/2Tq, where the threshold amplitude Tq is deflned by 
Eq. (25) in Ref. [^]. The above equations lead to the ratio /o//i = 3 and to the ratio of 
singlet to triplet cross sections 

^ = 9. (29) 

The beam target asymmetries (||) read CsTpp = 1, and Cbtpu = —0.8. The ratio of the 
total cross sections in pn and pp reactions is 5. Accordingly, the prediction for the spin 
density [|l^] reads poo = 0, pn = 0.5. 
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Let us now turn back to the Figs. ^ and ^. These figures show (i) a relatively small contri- 
bution of the direct radiation channel (b), which is in agreement with previous works P,p!5|, 
(ii) the cross sections for the pn interaction are qualitatively very similar in shape to these of 
the pp interaction but they are larger, and (iii) the ratio of the corresponding cross sections 
in pn and pp reactions is different for the sets A (or B) and C. Later we will discuss this 
aspect in detail. 

The energy dependence of the total cross sections of pp ppcj) and pn pncf) reactions 



for the sets B and C is shown in Figs. |I2| and [l3|. We do not display the result for the set 



A because it is practically the same as for the set B. The experimental data is taken from 



Ref. [|TT[]. One can see that the direct radiation channel (b) is much indeed smaller than 
the meson exchange contribution (a) in the near-threshold region where our consideration 
is valid. 

Fig. |l^ shows the energy dependence of the ratio of the total cross sections of pp ppcf) 
and pn pncj) reactions for the different parameter sets. One can see that this ratio increases 
with the energy excess and differs from the threshold value 5 in case of absence of FSI. The 
FSI is greater in the triplet initial (or singlet final) states and reduces the contribution of the 
initial singlet state in the pn interaction. For the set C the ratio o^nl^pp is greater because of 
the relatively greater contribution of the initial triplet state in the meson exchange channel 
(a). Fig. |15| shows the energy dependence of the ratio of the singlet to triplet cross sections in 
pn interactions (cf. Eq. (|29|)). The left panel shows this ratio for the separate channels, while 
on the right panel one can see our prediction for this ratio for the different parameters sets. 
One can again see a strong deviation from the threshold prediction (p9|) without FSI and 



non-trivial non-monotonic dependence of these ratios with some maximum values around 



As2 ~ 20 MeV. 



Fig. |T6| shows beam target asymmetry (^) for the separate channels for pp and pn in- 
teractions. For the pp interaction it coincides with its threshold value C^t = 1 up to a 
relatively large energy excess. For the pn interaction the asymmetry is different for the 
different channels which reflects the different role of the singlet and triplet states in the 
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different amplitudes which are additionally modified by the FSI. 

The total asymmetry for the different parameter sets is shown in Fig. |1^. It is interesting 
that even for the sets A and B with small contribution of the direct radiation amplitude (b) 
the asymmetry for pn interaction strongly deviates from the threshold prediction (without 
FSI: = —0.8), displaying a minimum around As5 ~ 20 MeV. 

We do not display here our results for the spin density matrix elements because for the 
sets A, B and C we get almost the threshold values, i.e. poo = 0, pu = 0.5, which reflects 
the dominance of the meson exchange channel (a). 

V. SUMMARY 

We have analyzed the production in nN and A^A^ interactions in the near-threshold 
region using the conventional "non-strange" hadron dynamics, based on the amplitudes 
allowed by non-ideal u — (p mixing, that is meson conversion in a (pTcp vertex and direct 
emission from the nucleon legs by a direct (pNN coupling. Using the limited body of 
available experimental data of the total unpolarized reactions we have tried to reduce as 
much as possible the uncertainty of the model parameters. As a result we get two branches 
of solutions with either a relatively small or a relatively large contribution of the direct 
emission channel which is determined by the strength of the (pNN interaction. By making 
use of these solutions we have compared various parameter sets with different strengths of 
the direct (pNN interaction. 

Analyzing the np —>■ ncj) reaction we find a strong dependence of the various observables 
on the strength of the (pNN interaction. The study of the differential cross section and 
angular distributions of electrons and kaons in the e~^e~ and —* K'^K~ decays seems 
to be most promising in investigating the (pNN dynamics. Experimentally, this study might 
be performed with the pion beam at the HADES spectrometer in GSI/Darmstadt. 

Analyzing the A^A^ —>■ NN(f) reaction we find a large difference in pp and pn reactions due 
to the different role of the singlet and triplet nucleon spin states in the entrance channel and 
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strong final state interaction. We predict a non-monotonic energy dependence of the ratio 
(^pn/(^pp and of tfie beam target asymmetry for tlie pn interaction wliicli deviates strongly 
from the pure threshold prediction. 

Finally, we emphasize once more that the present investigation is completely based on 
the conventional meson-nucleon dynamics and, therefore, our predictions may be considered 
as a necessary background for forthcoming studies of the strangeness degrees of freedom in 
non-strange hadrons. Additionally we would like to mention that fixing the (pNN coupling 
is important for an access to the elastic (pN scattering cross section which determines the 
degree of thermalization and collective flow properties of the mesons in heavy-ion collisions 
at SIS energies. 
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APPENDIX A: FINAL STATE INTERACTION 

In this appendix we present the formulae for the FSI and corresponding correction factors. 
We use the general framework for the FSI enhancement factor based on the Jost function 
formalism. Important aspects of this framework are described in the monograph by Gillespie 
27| and some early original papers [^. With respect to the significance of this problem 



in studying various near-threshold particle production reactions in the present time with 
cooled beams, we accumulate here the relevant expressions of this method and give the final 
result in a form convenient for specific calculations. For the Jost function formalism we use 
the notation of the textbook by Newton For simplicity, we limit our consideration to 
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the s-wave interaction which is dominant in the near-threshold region. A generahzation for 
higher angular momenta may be done straightforwardly. 

The FSI enhancement factor for two identical particles with momentum k in their c.m.s. 
reads 

where J'+ belongs to a set of functions J±{k) which are defined through the Wronskian of 
two linearly independent solutions of the Schrodinger equation, 

J±{k) = U{k, r) ^'{k, r) - r^{k, t) ^{k, r), (A2) 

where the prime means here the derivative with respect to r. The function J'= Jj^ is called 
the Jost function. The integral equations for the regular and irregular functions (f{k,r) and 
f±{k,r) have the standard form 

Sin AjT* 1 /'C'O 

(p{k,r) = — \- — / dr' sin k{r — r')V{r')f±{k,r'), 

k k Jq 

1 f°° 

Uk, r) = e^^'' - - / dr' sin k(r - r')V(r')Uk, r'). (A3) 
k Jo 

Eq. ( |A2D and the boundary conditions ^^(O) = 0, V''(0) = 1 show that J±{k) = f+{k, 0), thus 
allowing the integral representation 



J(k) = 1 + - / dr sin krV(r)f+(k,r). (A4) 
k Jo 

The physical wave function ip~^{k, r) is related to the regular function (p{k, r) as 

r(k,r) = ^^, (A5) 

which means that the inverse of the square of the modulus of J^{k) measures the probability 
of finding the particles near r = 0, relative to a situation without interaction. From Eq. ( |A4D 
one can find the important asymptotic condition 

lim J = 1, (A6) 

\k\—*oo 
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which shows that at high energies the enhancement tends to unity, thus leaving the total 
amplitude unchanged. 

The analyticity of J together with the asymptotic condition Eq. ( |X6| ) leads to the integral 
representation of the Jost function in terms of the phase shift 



dk 



(AT) 



k — k' + ie 

where k„ is related to the corresponding binding energy = — 2/ie„ > if bound states 
appear; is the corresponding reduced mass. 

For the practical usage of the above formalism it is convenient to work with the effective 
potentials which give an exact analytical expression for the phase shift. Let us first consider 
the singlet NN scattering (without bound state). The Eckart potential 



V{r) 



+ 



a^-(3^ \a-(3 a + (3 
gives the s-wave phase shift 



k cot oq = + 



a > 0,/3 > 0, 



a — (3 a — P 
reproducing the effective-range phase shift exactly, 

1 1 



k cot 6o 



ao 2 



rok 



with 



1 



1 



" = - ( V 1 - 2^o«o + 1 1 , /3 = - ( V 1 - 2roao ^ - 1 



The insertion of 6n from 



6o = -In 
" 2 



{k-ia){k + ip) 



{k + ia){k-i(3)_ 
in Eq. (|A7|) gives the analytical expression for the Jost function 

k + if3 



J{k) 



k + ia' 



(A8) 



(A9) 



(AlO) 



(All) 



(A12) 



(A13) 
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and the resulting enhancement factor for singlet interaction reads therefore as 

2x sin 60 e^^° 



Uk) = Coik' 



k 



Co{k' 



f.2 ^ ^2 {krof + 2(1- roflo ^ + v^l - 2roao ^ 



a — (3 



2rn 



(A14) 



which coincides with the classical Watson enhancement factor 

sin ^oe**^" 



w 



(A15) 



in the limit of kro 0. Note that expression (|A15|) is commonly used in calculations of the 
FSI in the near-threshold region (cf. |31|), where the constant Cw is fixed by a comparison 
of calculation and experimental data. Eq. ( |A14| ) is superior to Eq. ( |A15| ) because it yields 
a definite value of Cw = Co{k'^ = 0) by making use of the independent phase shift data and 
satisfies at the same time the required asymptotic behavior according to Eq. 



lim Io{k) = 1, 

|A;|^oo 



(A16) 



contrary to Eq. ( |A15| ), where Iw{k) at \k\ —>■ 00. 

Fig. |T8| illustrates the effect of FSI when using different enhancement factors for our 
parameter set B for the reaction pp pp4', where we keep here only the dominant mesonic 
exchange diagram in Fig. 2a and use the threshold value. In this case the energy dependence 
comes only from the enhancement factor and the phase space volume. The cross section 
calculated with enhancement factors from Eqs. ( |A14| ) and (|A15|) with Cw = Co(0) are 
shown by the dashed and the solid lines, respectively. For comparison we show also the 
cross section calculated without FSI (dot-dashed line). One can see, that the difference 
between the two factors ( |A14|) and ( |A15|) is indeed negligible at sufficiently small energy 
excess, say at As^ < 10 MeV with kr-Q <^ 1, where the Watson theory [BO] is valid, thus 



supporting the approach of |3T|. At As 2 > 50 MeV the difference between the two variants 
results in a factor 2 and greater. 

Following Ref. for the triplet pn interaction one can use the effective potential of 
the Bargmann type 
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V{r) 



-4k- 



d 
dr 



sh air 



g{K + ai,r) - g{K- ai,r)_ 



g{q,r) = (e '^^ + 2shqr)k 



(A17) 



which reproduces the known phase shift and the deuteron binding energy. The Jost function 
in this case reads 



J{k) 



k — IK 
k + iai 



(A18) 



with 



«! = (2 - Kri)r^ \ 



(A19) 



where = 2/ierf, is the deuteron binding energy, and ri and ai = — are the 

triplet effective radius and scattering length, respectively. The triplet enhancement factor 
reads explicitly 



Ii{k) = Ci{k') 



2. sin ^1 e*''^ 



Ci{k' 



p ^ ^2 {kn f + 2(1- na, 1 + - 2riar') 



ai + K 



2ri 



(A20) 



In our calculation we use ao,i and ro,i from Ref. 



pn singlet : aopn = — 23.768fm, rop„ = 2.75 fm, 
pp singlet : aopp = — 7.8098 fm, ropp = 2.767 fm, 

pra triplet : aip„ = 5.424 fm, rip„ = 1.759 fm, k"^ = 4.318 fm. (A21) 

Finally, we would like to mention that the approach presented here is equivalent to the 
approach in Ref. if the off-shell correction factor V (cf. Eq. (9) in |^) takes the form 
V{E, k) = -a^r-^ (l + ^1 - 2r^a-i) with x = 0, 1. 
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FIGURES 




FIG. 1. Diagrammatic representation of the ttN Ncj) reaction mechanisms: (a) meson 
exchange diagram with (p emission from the ^pTr vertex, (b) direct cf) emission from the (pNN 
vertex in Compton hke diagrams. 




FIG. 2. Diagrammatic representation of the NN — > NNc/) reaction mechanisms: (a) meson 
exchange diagram with cf) emission from the internal meson conversion in the (ppir vertex, (b) direct 
(p emission from nucleon legs. The zig-zag lines depict effective boson exchange. Exchange diagrams 
are not displayed. 
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FIG. 6. The same as in Fig. |5|, but for parameter set C. 
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FIG. 7. The differential cross sections for the ir p ^ ncj) reaction at As 2 = 50 MeV for the 
same parameter sets as in Fig. ^ 
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FIG. 8. The spin density matrix elements poo (left panel) and pu (right panel) for the different 
parameter sets at Asa = 50 MeV. 
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FIG. 11. The longitudinal (Wl) and transversal (Wt) fluxes of the outgoing electrons (left 
panel) or kaons (right panel) for various values of the anisotropies. 
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FIG. 13. The same as in Fig. |12|, but for the parameter set C. 
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FIG. 15. The energy dependence of the ratio of the initial singlet to triplet cross sections in 
■pn interaction for separate channels (left panel) and for different parameter sets (right panel). 



PP ->PP^ 



o 



1.2 



1.0 



0.8 



0.6 



pn ->pn(^ 





























\ \ 








\\ - 




\ N 








\\ 




\ \ 








w 




\ \ 



o 



10 

As'^' [MeV] 



100 



0.4 



0.0 



-0.4 



-0.8 



\ 

\ 

\ 

\ 

\ 

\ 


— - (a) 
— (b) 


- 

\ 

\ 




/- 

/ 



10 

As'^' [MeV] 



100 



FIG. 16. The beam target asymmetry for separate channels for (left panel) and pn (right 
panel) interactions. 
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FIG. 17. The same as in Fig. 15 but for different parameters sets 
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FIG. 18. The effect of FSI for different enhancement factors for the parameter set B in the 
reaction pp pp(j). Notation is explained in the text. 
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